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We study the reflection of an acoustic plane wave from a steadily sliding planar interface with 
velocity strengthening friction or a shear band in a confined granular medium. The corresponding 
acoustic impedance is utterly different from that of the static interface. In particular, the system 
being open, the energy of an in-plane polarized wave is no longer conserved, the work of the external 
pulling force being partitioned between frictional dissipation and gain (of either sign) of coherent 
acoustic energy Large values of the friction coefficient favor energy gain, while velocity strengthening 
tends to suppress it. An interface with infinite elastic contrast (one rigid medium) and V-independent 

■ (Coulomb) friction exhibits spontaneous acoustic emission, as already shown by M. Nosonovsky and 
G.G. Adams (Int. J. Ing. Sci., 39, 1257 (2001)). But this pathology is cured by any finite elastic 

f*"^ . contrast, or by a moderately large V-strengthening of friction. 

C^l ' We show that (i) positive gain should be observable for rough-on-flat multicontact interfaces (ii) 

(— i , a sliding shear band in a granular medium should give rise to sizeable reflection, which opens a 

■ promising possibility for the detection of shear localization. 



in 



45.70.-n, 43.40.+S, 46.55.+d 



. I. INTRODUCTION 



The question of the origin and nature of shear localization in disordered systems, such as soft glasses or confined 
granular media, which are jammed at equilibrium, but flow when sheared beyond a threshold stress, is a long standing 
one. Due in particular to recent progress in theoretical descriptions, it is presently the subject of renewed interest. 
£f) ' Hence the need for identifying appropriate, non invasive methods of experimental investigation which could comple- 
ment the optical and NMR imaging ones, recently put to use by Pignon et al [2], and by Coussot et al [3]. We intend 
to show in this article that the propagation of sound pulses appears as a promising possibility. Namely, we will show 
that the presence of a shear band of thickness small compared with the acoustic wavelength should give rise to strong 
anomalous reflection of a transverse acoustic signal for well defined incidence ranges. Such a method could therefore 
provide a relatively handy fingerprint of shear localization in confined granular media. 

An extreme case of localized shear flow is that of frictional sliding of the interface between two cohesive macroscopic 
solids. In such a configuration, the very structure of the system prelocalizes shear to the nanometer-thick layer which 
forms the molecular adhesive contact (s). The role of the above-mentioned threshold stress is played here by the so- 
called static friction force. Below this threshold, the interface is elastically pinned (jammed), and responds elastically 
• i-H , to a shearing force. So, the corresponding mechanical boundary condition is simply that the displacement fields in 
the two media must be fully continuous across the interface. But, beyond this shear level, sliding sets in and, along 
the direction of relative motion, the boundary condition is now provided by the dynamic friction law, which states 
that the shear and normal stresses are proportional. Obviously, such a discontinuity of the boundary conditions must 
result in a discontinuous change of the acoustic reflection and transmission coefficients when the system is set into 
sliding. 

This was pointed out already long ago by Chez et al [4] , who studied the reflection of a sound wave propagating 
in a plane orthogonal to the sliding direction, and polarized in the plane of incidence on an interface with the 
Coulomb frictional behavior (constant friction coefficient). However, due to the choice of this particular geometry, 
they overlooked an interesting effect. Indeed, the sliding system is an open one : energy is being pumped in from 
an external source - the external driving machine which imposes the sliding velocity. So, as soon as the acoustic 
displacement field has a non-zero component, in the interfacial plane, along the direction of the sliding motion, 
additional mechanical work (of a priori either sign) is extracted from the external source, and interfacial acoustic 
scattering is no longer energy-conserving. This opens in principle the possibility of acoustic gain at reflection, i.e. 
conversion of incoherent into coherent mechanical energy - quite an exciting prospect indeed! 

Now, from the point of view of the propagation of an acoustic signal of wavelength A, a shear band of thickness 
d << A in a confined granular medium appears as the equivalent of a frictional interface between two identical 
solids. Indeed, the band can then be considered as a surface of mechanical discontinuity between the non sliding 
adjacent regions, which behave as (disordered) elastic solids. Experimental studies by rock mechanicians [5] [6] of 
systems constituted of two bulk rock pieces separated by an interposed layer of granular material (called "gouge") 



1 



have established that in such systems (i) sliding occurs in a narrow band within the gouge (ii) the dynamics is ruled 
by a standard solid friction law, the associated friction coefficients having a magnitude comparable with those for 
solids in direct contact. 

Reflection and transmission of waves with a polarization component along the sliding direction of an interface 
with constant friction coefficient between dissimilar media have been recently studied by Nosonovsky and Adams 
[7], though in a different perspective. Namely, they focussed primarily on the possible generation of slip pulses - a 
dynamical feature that seems to be specific of pure Coulombic friction. In the course of this article we will rederive 
some of their results, which will be extended to the more realistic case of velocity-dependent friction and to the shear 
band problem.. 

This article is organized as follows. 

In Section 2, we first write down the equations for the most general case of a monochromatic acoustic wave incident 
upon a planar frictional interface between two semi-infinite solids with different elastic moduli. We then specialize 
in Section 3 to the case where the elastic contrast between the two media is very large (e.g. a gel sliding on top 
of a glass plate). We show that, if the stiffer medium is assumed strictly non dcformable and friction taken to be 
Coulomb-like, the reflection coefficient of a wave with the sliding direction and polarization in the plane of incidence is 
highly pathological : not only is a huge gain at reflection possible for some particular incidence range, but spontaneous 
acoustic emission from the surface is predicted - a result already obtained recently by Nosonovsky and Adams [8]. 
These singularities are to be related with those already found by Adams [9], and Ranjith and Rice [10], in their 
studies of interfacial waves in the same system. They result, as is well known in mechanics, from the specific singular 
character which the Coulomb model, which takes the friction coefficient to be a mere constant, shares with the Hill 
model of plasticity. Indeed, we show that (i) a very small finite relative compliance of the stiffer medium is sufficient 
to destroy the acoustic emission singularity (ii) improving upon the Coulomb description by taking into account a 
velocity-strengthening dependence of the dynamic friction coefficient of the order of what is measured for real systems 
also cures this singularity. Moreover, the possibility of energy gain at reflection is found to be strongly dependent 
upon the strength of the velocity dependence, and hence on the type of system : while it should be negligible for 
a gel-on-glass system, it might be observable with a rough-on-flat multicontact interface in well defined incidence 
ranges. 

Section 4 is devoted to the symmetric case of two mechanically identical solids relevant to the shear band problem. 
In this situation, and when the velocity dependence of the friction coefficient is taken into account, we predict that 
the sliding shear band should give rise to a clear acoustic signature - namely a reflection coefficient with magnitude 
of order typically of order unity for incident signals in well defined ranges of incidence angles. 



II. GENERAL FORMULATION 



Consider two elastically isotropic semi-infinite media (M; M') with shear moduli and Poisson ratios (fi, v ; p! , v') 
and densities (p, p'), occupying respectively (Figure 1) the upper {x 2 > 0) and lower (x2 < 0) half spaces, and in 
continuous contact with each other along the x 2 = plane. Medium (M') is assumed to be in stationary sliding 
motion with respect to (M) at velocity v along X\ and towards x\ < 0. In this reference state, the (homogeneous) 
normal and shear stresses t 22 , t± 2 , are related by the dynamic friction law : 

r*, 2 = -f{v)T* 22 (1) 

with f(v) the friction coefficient. 

An emitter linked to (M) is sending from infinity towards the interface a plane acoustic wave of frequency u, 
propagating in the xi,x 2 plane at incidence angle 9 (Figure 1) and polarized in the plane of incidence. That is, the 
associated displacement has a finite component along the sliding direction. In order to fix ideas, and for the sake of 
simplicity, we restrict in most of what follows the algebraic formulation to the case of a transverse (shear) incident 
wave - the case of an incident longitudinal (dilatational) signal follows straightforwardly. The elastic displacement 
field (ui,u 2 ) in (M) obeys the Lame equations : 

o d 2 u\ , 9 9 , d 2 u 2 , d 2 U\ , . 

U ^ = C ^ + ^-^3^t 2 +C -34 < 2 > 

2 d2u 2 , 2 2 \ <9V 2 d 2 U 2 

U ^ C ^ 2+{CL - CT) dx^dx- 2 +CT dxV 2 (3) 
cl,ct are the longitudinal and transverse sound velocities in (M), and 
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4 = ^ % = ^^ (4) 

The solution of equations (2), (3) reads : 



u 2 = eeUkxi-wt) 



k 



q T q T k 

where e specifies the amplitude of the incident wave, k is the x\ component of its wavevector, and 



(6) 



1T.L - j 



k 2 



(7) 



T,L 



with qr real positive and Req^ > and Im^L > 0. 

In the lower medium (M'), which moves at velocity (—v) in our reference frame, analogously : 

u[ = ee i( - kx i-"V \ a 'e- iq T X2 + e~ iq! ^ 



(8) 



/ = ee i(k Xl -u,t) 



q T k 



(9) 



where 



/2 

J = cj + vk q' TL 2 = k 2 (10) 

c 't,l 

In order to determine the amplitude reflection (a, (3) and transmission (a' , /?') coefficients into the transverse and 
longitudinal channels we must now specify the four boundary conditions (BC) to be satisfied along the deformed 
M/M' interface. 

Since the acoustic stresses we consider are small, 0(e), we expect the contact to persist everywhere, so that 

(i) normal displacements on both sides of the interface are equal : uh — u'h and mechanical equilibrium imposes 
continuity of normal and shear stresses : 

(ii) i~ nn t nn 

(iii) T n t = T' n t 

The sliding velocity perturbation (vi — v) being also 0(e), vi remains positive (sliding persists) everywhere, and 
we assume that the dynamic friction condition now holds locally, namely that, at each intcrfacial point : 

(iv) T nt + f(vi)r nn = 

These boundary conditions, which have been taken for granted in existing studies of interfacial waves and shear 
fractures, actually imply important physical assumptions, which should be made explicit. 

On the one hand, the (macroscopic) contact is supposed to be continuous and homogeneous - hence the statement of 
homogeneous stresses in the reference state. Now, it is well known that such is not the prevalent case. Most interfaces, 
being formed between rough solids, are actually constituted of a large, but rather sparse, set of microcontacts [11] 
[12], in between which the surfaces are mechanically free. Then, it is clear that conditions (i)-(iii) are valid only in 
a coarse-grained sense, i.e. provided that the length scale 2n/k of variation of the acoustic fields along x\ is much 
larger than the average distance d between microcontacts, that is for fee? << 1. 

Intercontact distances lie commonly in the range of hundreds of micrometers. Only in the case where one at least 
of the solids is extremely compliant are macroscopic contacts truly continuous. This is for example the case for 
elastomers or gels. 

On the other hand, the mere existence of a finite static threshold proves that frictional dissipation results from the 
triggering by the applied shear of structural instabilities. It is now documented that the corresponding structural 
rearrangement events take place in the nanometer-thick adhesive interfacial layers, and affect volumes of, typically, 
nanometric scale [12] [13] - comparable with the "shear transformation zones " invoked by Falk and Langer [14] to 
model the plasticity of amorphous solids. Let us call b the size of such a zone. A friction law represents the result of 
a statistical average over a large number of such dissipative events. So, it can only make sense on a scale much larger 
than b, that is for kb « 1. 
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All this means that boundary conditions (i)-(iv) above must be understood as valid only on scales larger than some 
cut-off length L — max(6, d) on the order of, typically, (i) a fraction of micrometer for conforming solid contacts (ii) 
a fraction of millimeter for the, more common, multicontact interfaces. 

Finally, since Vi now becomes a time-dependent quantity, condition (iv) implies that we assume the friction relation 
to hold instantaneously on the scale uu^ 1 . It is known that this does not apply to the case where the steady state 
f(v) is velocity weakening (df /dv < 0). Indeed, such behavior necessarily results from the action of some underlying 
structural dynamics leading to aging when sticking and rejuvenation upon sliding, such as that associated with the slow 
creep growth of the real area of contact relevant to multicontact interfaces [15] [16]. Then, the equations describing 
non steady friction must involve explicitly at least one more dynamical "state" variable, and condition (iv) becomes 
insufficient. Note also that it is in this regime that steady sliding may be unstable with respect to stick-slip. 

So, we restrict ourselves in what follows to the velocity-strengthening case df /dv > 0. This can be expected to 
hold, for rough-on-rough interfaces, only at sliding velocities in the mrn.sec -1 range, large enough for rejuvenation 
effects to be saturated [17]. However, it has been shown [18] to prevail down to the /xm.sec -1 range when working 
with rough-on-flat interfaces where contacts keep their identity when sliding, which makes contact area saturation 
easily realizable. 

The position of our deformed interface is given by 



X 2 I (xi,t) = 1t 2 {xi - Ui (xi,0,t)) 



(11) 



We assume from now on that acoustic deformations are small enough for us to work in the linearized approximation. 
Then, the normal and tangent unit vectors n, t are simply : 



L,*)=^ ^■■^ r ,0, 



1 



«,(.ri../)= j V /. n .i)./ j t(xi,t) = y ^du 2 ^j 

Let the stress field in, say, M, be denoted t*j + STij, (i,j = 1,2), with 

Then, to first order, at the interface in M 



(12) 



dui duj 2v 

Q 1" "5 1" i 7^"ij u kk 

axj oxi 1 — 2v ■ 



TliTi-jTlj — T 



22 T "722 



St 22 - 2Ti*o 



( dU2 

\dxi 



(13) 



(14) 



and the local sliding velocity 



Conditions (i)-(iv) then become: 



T n t = T~i2 + ST12 + (T 2 * 2 - T * x ) 

vi = v + (iii - ii'i) 

[«2 - u' 2 ] xifi t = 

St 22 - St 22 - 2t* 2 



( du2 
\dxi 



( du 2 




\dxi 


dxi) 



= 



(15) 

(16) 

(17) 
(18) 



xi,0,t 



St 12 - St[ 2 + (t 22 



22 T ll) 



/ du 2 




\dxi 


dxi) 



(19) 



6r 12 + f(v)5r 22 f'(v) T * 22 (in - u[) 



(20) 



xi,0,t 
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where f'(v) = df/dv, and we have set 



(r 22 



2f(v)r-_ 



12) 



(21) 



r 22 < (compressive normal stress), Tj* 2 = — r 22 > , and we can reasonably assume sliding to occur under zero 
lateral stress, i.e. t*^_ — 0, so that, in general, r* > 0. 

Note that the r*-dependent terms in eqs (17)- (20), which account for the fact that the BC's must be enforced along 
the deformed interface, have usually been overlooked in previous works on interfacial waves. As will appear below, 
they give rise to corrections on the order of (r*//x). In the case of hard materials, externally applied stresses are 
in practice always considerably smaller than elastic moduli, and these corrections can safely be neglected. However, 
such may not be the case when dealing with very compliant materials such as gels or elastomers. For example, for 
the case of gel/glass interfaces [19] sliding stress levels may be a sizeable fraction of /x, and the full expressions should 
be retained. 

Then, using eqs (5), (6), (8), (9) and (13), the BC's yield the following set of four non homogeneous linear equations 

for a, a' , (3, f3' : 



k . Ql R k , q' L , k 

a+— p T a + — p = 

q T k q T k q T 



(22) 



-2/xa + /! 



Qt 



k 2 



— /3 + 2/A*' -» ,q ' T 



k 2 



-0 = 2/i 



(23) 



q% 



q T 



k \ + 2 mL p- + ^rt T ^a' + 2ti>qi L {3> = ^ e 



q T 



q T 



(24) 



Qt 



k 2 t* k 2 

q T fi q T 



2q L + I 



Q% 



+^(a + p-a' -0') 



Qt 



—Ql P 
A 4 J 

k 2 of , T* k 2 

— + 2fk + 

q T H qr 



(25) 



where we have defined : 



A = 



f(v) CT 



f(v) 



(1 



(26) 



which measures the dimensionless strength of the velocity-strengthening effect. 

Equations (22)-(25) above are relevant to the case of a transverse incident wave. We will also display some results 
for a longitudinal one. In this latter case, the only modification concerns the first terms in the r.h.s. of equations (5), 
(6) which give the expression of the total acoustic displacement field. This results in leaving the l.h.s.'s of the final 
equations (22)-(25) unchanged, only the r.h.s.'s being modified. 

Since the general solution of our problem depends on a large number of parameters (four elastic moduli, two 
densities, the angle of incidence, the friction coefficient and its derivative) it will be more illuminating to focus on a 
few simple cases, namely that of very large elastic contrast (quasi non deformable M'), and the symmetric case of two 
identical materials, relevant to the problem of shear band detection. 



III. LARGE ELASTIC CONTRAST CASE 



We consider here the case of maximum asymmetry, where M' is considerably stiffcr than M. In order to disentangle 
the respective effects of the elastic contrast ratio R = it'/ yU and of the velocity dependence of friction, we first treat the 
extra simple limit of a non deformable M' and of pure Coulomb friction, then examine how the results thus obtained 
are modified (i) when R is large but finite (ii) when f'(v) ^ 0. 
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A. Non deformable M' and Coulomb friction 



The boundary conditions reduce to imposing that ui = and T nt + fr nn = along the xi = plane. Equations 
(22)-(25) become, in this R- 1 = f = limit, where a' = /?' = : 



—a- —0 = — 

q T k q T 



(27) 



k 2 -q 2 



q T 



2fk 



-2q L + f- 



qr 



k 



qf 



q T 



2fk 



which yield for the T — > T and T — > L reflection cocfflcients : 



where : 



a = 1 + 



qL 



q T A 



q T k 



f[q T - ^qh 



q T 



(28) 



(29) 



(30) 



Note that, for / — > 0+, and for incident waves propagating towards X\ > (k > 0), A > 0. That is, at least in the 
limit of weak friction and in this incidence range, clearly a > 1, (3 > 0, and the reflected energy flux is larger than the 
incident one. In order to try and clear up this a priori surprising prediction, let us look in more detail into the question 
of energy balance in our system. Consider the volume (V) of a slice of unit depth along X3 of medium (M), limited by 
its area (S) {—L < X\ < L) along the interface, and the semicylinder (Coo) of radius L(L — ► 00). In the situation we 
are considering, of an incident wave of constant amplitude, the average over an acoustic period of the elastic energy 
stored within (V) is time independent. This means that energy conservation simply imposes that the net energy flux 
at infinity Woo, flowing out of Coo, associated with the acoustic waves, must equal the increment associated with the 
acoustic perturbation of the work injected per unit time into the system via the work of the stresses acting on the 
moving interface, W ext , which is pumped from the driving machine. This is proved in detail in the Appendix, where 
we show that, in the present case on a non deformable M', per unit area of (S): 



(Wext) = (-M1T12) 



(31) 



where (. . .) stands for the average over the acoustic period. We then easily prove (Appendix A) that this is exactly 
equal to the net acoustic energy flux : 



(Wrefl) ~ (Wi, 



(32) 



In other words, a sliding frictional interface is in principle able to transform external incoherent mechanical work 
into coherent acoustic radiation - i.e. to act as an "acoustic laser". In order to make this more precise, we plot on 
Figure 2a the relative powers reflected into the T and L channels, given by (see Appendix): 



wtt =\ a 



(33) 



WTL = 



qrqL 
k 2 
= 



(3 I if qL real 

if q L imaginary 



(34) 



against the incidence angle 9 = sin -1 (crk/uj) and for / = 0.2, ct/cl = 0.5. It is seen that, for all positive incidences, 
a transverse incident wave is predicted to be notably amplified upon reflection, while the power in the L channel 
remains quite small, even close to 9u m = siu^ 1 (ct/cl) beyond which the L wave becomes evanescent, where it is 
maximum. However a much more surprising result, already derived by Nosonovsky and Adams [8] is that wtt and 
wtl are found to diverge at a negative incidence — cri where the denominator A (Eq.(30)) of the expressions for a 
and (3 (Eq.(29) vanishes. Using equation (7), the condition A = reads : 



'^-sin 2 



/cos 29 sin 9 
1 - /sin 29 



(35) 
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which is easily checked graphically to have a single negative solution — 9 cr with 9 cr < 6u m , whatever the values of / 
and of ct/cl- The smaller /, the closer 9 cr approaches 9u m . In the case shown on Figure 2a, 9u m — 9 cr ~ 7'. In 
the cease (Figure 2b) of L incidence, the corresponding singularity occurs for 9 close above —90° In other words, in 
this admittedly oversimplified limit (rigid M', pure Coulomb friction), homogeneous sliding at constant velocity is 
impossible : indeed, any infinitesimal perturbation is able to trigger the emission, all along the interface, of a coupled 
set of transverse and longitudinal acoustic waves, of a priori undetermined amplitude, travelling respectively at angles 

9 cr and sin -1 c 1 sin^ e ^ m tnc back direction — x\. On the one hand this entails that an infinite external 

energy would be pumped in. On the other hand, as soon as the local interfacial velocity will vanish, the interface will 
stop (stick) . So this pathologic behavior signals that homogeneous sliding is here absolutely unstable - a signature to 
be added to that provided by the existence in such systems of amplified interfacial waves, first identified and studied 
by Adams [9], which lead to the question of ill-poscdness of the problem of interfacial slip pulses studied by several 
authors and recently synthetized by Ranjith and Rice [10]. Whether or not the stable mode of sliding for our model 
interface would be a set of square slip pulses such as calculated in ref. [8] is out of the reach of this work. We will 
rather concentrate on a different question, namely : how robust is this pathology? Does it persist in the presence of (i) 
a finite elastic compliance of the stiffer medium M'; (ii) a velocity-strengthening dependence of the friction coefficient? 
We now consider successively these two possibilities. 



B. Finite elastic contrast and Coulomb friction 



Let us now consider the case of a large but finite elastic contrast ratio R = fi'/fj, » 1. We assume for the sake 
of simplicity equality of the Poisson ratios v = v' . In this case, except for quasi- normal incidence angles, q% L < 
: there is total reflection at the interface, with only evanescent transmission into M'. a and (3 must now be obtained 
by solving the full equations (22) - (25) for constant /. Except in the vicinity of 9 = — 9 cr , a finite R^ 1 simply acts 
as a regular perturbation, leading to a correction O For 9 « — 9 cr , however, one must consider in detail the 

behavior of the determinant V of the system (22) - (25). To first order in R^ 1 : 

V cx A (q T , q L ) + R^T (q T ,q L , q' T , q' L ) (36) 

As seen above, the zero of A always occurs for 9 cr < 9u m where q^, and hence A, is real. On the other hand it 
appears that, for q' T L pure imaginary, T is a complex quantity. So, no zero of V which would evolve continuously 
from — 9 cr exists, and one finds that, for finite R -1 , a and (3 no longer diverge, but exhibit, in the vicinity of — 6 cr , 
a maximum of order R : a finite compliance of M', however small it is, is sufficient to kill the pathology found in the 
totally rigid limit. This we have confirmed by computing the relative reflected powers wu (I, J = T, L) for various 
values of R. The results for wtt are shown on Figure 3: the smaller R, the lower the maximum of w, which reduces 
to a few units for R < 40. Of course, in practice, for R >> 1, the very large values of the reflection coefficients at 
maximum mean that a very small incident amplitude e would suffice to induce a finite response, such that v + iii 
would vanish, leading to interfacial stick - i.e. homogeneous slip, though no longer absolutely unstable from a strict 
mathematical point of view, would be very weakly stable, as it could not resist perturbations of finite but very small 
amplitude. 



C. Infinite elastic contrast and ^-strengthening friction 

The problem is now specified by eqs. (22) and (24) in which a' = (3' = 0, and the position of the singularity of a, (3, 
if any, is therefore given by the zero of: 

A A = q L (^ + f) + f(q T -2q L -^) + A^M(f + ±-) (37) 
\q T k J \ q T J c T \k q T J 

where the dimensionless parameter measuring the velocity dependence effect is defined by equation (26): 



A = f'(v)cT I t 2 *2 I 
f(v) V 

When solving numerically for A^ = 0, we find that the singularity disappears for A > A cr . More precisely, we find 
that, as A increases, the zero of A^ approaches —9n m , which it reaches for A = A cr , beyond which it disappears, 
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due to the fact that no zero can occur in the regime where the L wave is evanescent. The threshold value A cr is 
independent of the value of f(v), but it depends upon the sound velocity ratio. For a reasonable choice of this 
parameter (ct/cl < 0.7), A cr is at most a few units. For example, for the presently used value ct/cl — 0.5, we 
obtain A cr = 1. This behavior is illustrated on Figure 4, where we plot the relative reflected power wtt against 
incidence angle for various values of A. It is also seen that the larger A, the smaller the gain at reflection. It appears 
that at another characteristic A value (A = 2 for our parameter values), the gain becomes negative (wtt < 1) for all 
incidences. That is, a strenghtening v-dependence of / is very efficient to kill the amplification effect. 

It is therefore important to evaluate an order of magnitude of A for real interfaces with a large elastic contrast, a 
good example of which is that of gel/glass couples. Frictional sliding at the interface between glass and a 5% gelatin 
aqueous gel has recently been studied in detail by Baumberger et al [19]. Using their data, we find that, for v in 
the mm.sec -1 range, f'(v)/f(v) w 2.10 3 seem -1 , while (c T \ t% 2 \ / V) ~ 2 m.sec -1 , so that A w 4.10 3 is a very large 
number. 

In this A >> 1 limit, the rigid M' version of equation (24) reduces, to lowest order in A -1 , to : 

a + 13 + 1 = (38) 

that is, due to the high relative cost of increasing its instantaneous velocity, the interface remains locked to the 
homogeneous sliding state and u\ = 0. 

So, for gel/glass sytems, acoustic reflection should not in practice be able to distinguish between the static and the 
sliding interfaces. 

Another realization of the high elastic contrast situation is provided by the multicontact interface between rough 
glassy PMMA and atomically flat silanized glass, recently studied by Bureau et al [18]. For this couple, p! j [i — 20. 
Thanks to the flatness of the glass surface, it is possible to saturate the slow growth of the real area of contact, which 
is responsible, in the case of rough-on- rough systems, for the velocity weakening behavior of the dynamic friction 
coefficient and of the associated stick-slip dynamics. Then, for v > l^m.scc -1 , f(v) is velocity strenthening and of 
the form : 



f(v) = fo 



1 + Clog 



(39) 



where / ~ 0.2, and £ « 2 — 4.10 2 . With Upmma = 1 GPa, ct — 10 3 m.sec 1 , and under normal stresses on the 
order of 5 kPa, one gets : 

A. 100 ' 200 (40) 

u fim.sec 1 

That is, for sliding velocities in the 100 ^m.sec -1 range, A is typically on the order of 1 — 2. 

Finally, note that, since in this configuration the average distance between the micrometric regions which form the 
real area of contact is on the order of a fraction of millimeter, such an experiment would ask for acoustic signals in 
the range of of a few hundred kHz at most, in order for our continuum description to be valid. 

The relative powers reflected into the two channels for T and L incident waves under these conditions are plotted 
on Figures 5 and 6, for A = 0, 1, 2. (Since in this case typical values of t* / fi and v/ct are < 10 -5 , calculations have 
been performed for r* = 0, w' = w.) Inspection of these results and comparison with the case of the non moving 
interface (Figure 7) lead us to the following conclusions. 

• The more favorable channels for observing gain at reflection are the diagonal (LL and TT) ones. 

• In both cases, a very narrow peak of w, reminiscent of the singularity obtained in the limit of Coulomb friction 
on a rigid substrate, is predicted for 9 ~ —6u m . However, due to its narrowness, observing it would be certainly 
be excessively demanding in terms of directional accuracy. 

• It thus appears more feasible to investigate one of the two following configurations : LL at large negative 
incidence angles, in the range of —45°, and (TT) at large positive incidence 9 « 60°. 

• Gain decreases very rapidly with increasing A, i.e. with decreasing v (see eq.(40)) : while, for A = 1 it reaches 
up to 20% in (LL) and 35% in (TT), already for A = 2 it has practically collapsed to zero. So, according to 
equation (40), the best situation corresponds to the largest possible driving velocities, in practice in the range 
of a fraction of mm/sec. 

With these conditions in mind, it seems quite possible to obtain substantial acoustic gain at reflection on a sliding 
rough-on-flat multicontact interface. 
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IV. SYMMETRIC CASE 



We now turn to the opposite limit where M and M' have the same elastic properties. It is a priori relevant to the 
case of multicontact interfaces between two pieces of the same material. However, the above mentioned rough-on-flat 
configuration is inappropriate in this case, since the asperities on the rough surface would give rise on the flat one 
to unavoidable and poorly qualified plastic damage and wear effects. On the other hand, as already mentioned, in 
the rough-on-rough configuration where such problems are irrelevant, geometric aging interrupted by motion results 
in a velocity-weakening behavior of f(v) for velocities up to at least the mm. sec -1 range, not very easy to acess in a 
stationary sliding laboratory experiment. So, such symmetric solid on solid interfaces turn out not to be well adapted 
in practice to test our predictions. 

The most interesting case, as already mentioned in Section I, is that of an established shear band, in particular in 
a highly confined granular medium. In such highly disordered systems, internal stress inhomogeneities (the so-called 
stress-chain phenomenon) of course give rise to scattering of acoustic waves. However, this is all the smaller that 
the effective acoustic wavelength in the medium is larger with respect to the correlation length of stress fluctuations, 
known to be on the order of a few grain diameters D only. In this kD << 1 limit, as shown by experiments on 
the propagation of acoustic pulses [20], where one can separate out unambiguously the signal corresponding to the 
propagation of a "coherent pulse", a description in terms of an effective acoustic medium is legitimate. Note that, as 
shear band thickness is also expected to be, for low shearing rates, on the order of a few D, kD << 1 also ensures 
that we can safely neglect it when dealing with our reflection-transmission problem. 

Laboratory experiments on the frictional behavior of layers of granular rock confined between granite plates [5] 
[6], and of glass bead assemblies [21] have shown that sheared confined granular media obey standard friction laws - 
namely, beyond a static threshold, T12 = —fr-xi- However, the velocity dependence of the dynamic friction coefficient 
is not yet very clearly established. The data on rock point towards weakening at low velocities. However, they are 
strongly dependent upon the level of humidity. This, together with the very small size of the grains used in these 
experiments, strongly suggests that capillary condensation around the intergrain Hertz contacts is responsible for a 
slow strengthening of the medium, interrupted by sliding, which should become negligible either at low humidity or 
for larger grains. On the other hand, in the absence of such slow transients, since granular systems are essentially 
athermal, one expects logarithmic dynamic strengthening of the type found for multicontact interfaces (see equation 
(40)), which results from thermal activation effects [12] [13], to be completely negligible. If such is the case, / would 
be v- independent in the low velocity, quasi static sliding regime. Note that this is what has be found by Geminard 
et al [21] for glass beads completely immersed in water - though under relatively weak confinement. 

So, let us first consider the simplified case of pure Coulomb friction (A = 0). We then solve the corresponding version 
of eqs.(22) - (25) where we neglect t* /p, since, for the confined granular systems made of hard materials (such as glass 
or steel) which we have in mind, under ordinary experimental conditions [20] [22], this does not exceed, typically, 10~ 5 
at most. Moreover, as we are interested here in strongly subsonic sliding velocities (v/c T < 10~ 5 — 10~ 6 ), ui' ~ u. 

The results for the four relative reflected powers are plotted on Figure 8 for the case / = 0.2 : while, in the (TT) 
channel, for negative 9 there is loss at reflection, wtt is seen to be sizeable everywhe re in the range ( 0iim]®iirn) 
where a propagating L reflected wave exists, as well as for quasi-grazing conditions. For 8 = ±8u m the reflectance 
curve exhibits cusp-like maxima - a standard behavior for multichannel scattering cross sections at the closing point 
of a channel. About 6 = 0u m , we even predict a gain of more than 50%, which becomes much larger for / = 0.4. In 
the other three channels, w is much smaller - except, for the LL one, in quasi-grazing conditions which are hardly 
realizable in practice. 

This behavior is to be contrasted with what is expected from non localized homogeneous shear sliding, namely 
undisturbed propagation all through the sample. So, one expects the following to happen after starting to shear a 
confined sample at a constant rate, if a transverse acoustic pulse is sent into the granular pack in the direction normal 
to the shearing plane. At very short times, when the stress has not yet reached the sliding threshold, the system 
experiences uniform elastic deformation, the incident pulse will give rise to a "coherent" reflected one, which will have 
travelled way and back across the whole container. After the initial sliding transient corresponding to the gradual 
installation of the shear band, during which fluctuations are expected to kill the coherent signal, the reflected pulse 
should reappear, but with a distinctly shorter travelling time, thus providing a clear signature of shear localization. 
Note that, since its travelling length is reduced, so will be the attenuation due to scattering by stress fluctuations, 
making it more easy to detect than the former signal. Finally, when sliding is stopped, the situation should revert 
practically to that before shearing, though the density contrast Sp/p, on the order of a few % at most, associated 
with the shear band [23] still persists. Indeed, then, wtt ~ (Sp/ p) 2 (oIuj/ct) 2 is negligibly small, and pulse reflection 
only occurs, again, at the bottom end of the container. 

The question is then to check how robust the reflectance characteristics predicted for A = are with respect to 
a possible, though probably small, velocity strengthening dependence of the friction coefficient. It is seen on Figure 
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9 that, although, as expected, wtt decreases as A grows, at normal incidence it remains non negligible up to the 
sizeable value A — 2, where it still is on the order of 25%. 

V. CONCLUSION 

In summary, from the above results, we conclude to the interest of performing experimental studies of the reflection 
of acoustic pulses on sliding solid interfaces and shear bands, the potential interest of which is different in each case. 

As discussed in section III.C, experiments with a rough-on-flat multicontact interface open the possibility of obtain- 
ing acoustic gain, through the conversion by the frictional system of incoherent mechanical energy pumped from the 
driving system into coherent acoustic vibrations. However, such experiments are certainly delicate to realize, since 
they ask for working at non normal incidences, and, due to the effect of velocity strengthening of dynamic friction 
coefficients, at not very small velocities, in the range of typically a few hundred /xm/sec. 

On the other hand, acoustic pulse reflection appears as a promising method for detecting shear localization in a 
confined granular medium. The best configuration - a transverse pulse at normal incidence - is more easily realizable 
since it may be implemented with a single transducer acting as both emitter and receiver. Moreover, the expected 
signature (a strongly reduced transit time before return of the reflected pulse) should be easy to identify - its presence 
confirming at the same time the expected quasi-independence of the sliding stress on the sliding velocity. 

We therefore strongly hope that this work will motivate such experiments in the near future. 
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APPENDIX A: ENERGY BALANCE AT A FRICTIONAL INTERFACE 

We restrict ourselves to the case of a rigid substrate with an infinite elastic contrast R = ■— > oo considered in 
Sec. Ill A. This situation is sketched in Fig. 10, which concerns specifically the case of a transverse wave incident 
below the critical angle 6u m , so that also the longitudinal reflected wave is propagating. In our 2D geometry, the 
incident beam irradiates a strip So of the interface perpendicular to the sliding direction x\ and extending from — Lq 
to L along this direction. It has an area So P er unit length along x 3 . The reflected beams stem from the irradiated 
area. The beams are wide enough to make the fringe effects negligible and easy to exclude as additive constants to 
the principal quantities proportional to So- The So strip is now overlapped by a wider region S stretching from — L 
to L such that a semicylindrical dome Coo raised over it encloses the regions of beam overlap and interference (as well 
as the fields of attenuated waves, should they arise). The bottom of the dome is infinitesimally above the interface, 
so that the equations of motion (2), (3) are valid both inside and on the surface of the region V enclosed by the dome 
and the bottom plane. Then, the total elastic energy W inside the dome can be studied using the differential energy 
conservation law [24], [25] (see also [4]) 

divj + e = i.e., je t t + e = 0, (Al) 

where e is the local energy density, whose form will not be needed, and the energy current density j is given by 

3 = -ut i.e., j e = -UjT je (A2) 

We consider elastic fields having a steady homogeneous component and a single frequency acoustic component. As 
usual, in bilinear expressions like Eq. (A2), we turn to real parts of all quantities involved, which is consistent in a 
linear theory. Thus, we have 

u(r, t) = Rc(-iuju°(r)c- luJt ) (A3) 

x(r, t) = T* + Rc(ST°(r)e- luJt ) (A4) 
Introducing this into (A2), we may single out the steady flow component of j by averaging over the period 2tt/lu: 

(j) = ±(-iuju°'5t 1 > + iuJu^dr ) (A5) 
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where overbar denotes complex conjugation. The time average (W) of the total elastic energy per unit length of the V 
region for a stationary irradiation is zero. Integration of the time average of Eq. (Al) then yields a balance equation 
for the steady time averaged energy flows across the surface of V: 

J dS(j)-n = J dS(-(j»-n (A6) 

Coo s 

with n the outer normal of V. With the sign conventions of Figure 10, we obtain the following relation for the gain 
Woo, defined as the imbalance between the incident and reflected energy flows: 

Woo = {W T ) + (W L ) -(W inc ) = (W ext ) (A7) 
s ■ ' 

(Wrefl) 

The energy flows associated with the plane wave beams crossing the dome Coo are easily obtained, if we use the well 
known steady current density for a plane wave [24] 

\(j)l,t\ = ^ 2 cl,t-\U\ 2 

valid for both polarizations; the quantity U is a possibly complex amplitude. The plane waves appearing in (A7) 
are defined in Eqs. (5), (6). The coefficients a, (3 were obtained in Sec. Ill A. The three total energy flows are the 
products of current densities and beam cross sections. Thus, 

(w mc ) = >v- e 2 (— V-So^ 

\c T q T J uj 

<Wt) = 1„A T -«>P(;^ 2 -So^ (AS) 

The ratios of the flows wtj = (Wj) / (W n c) {J = T, L) then have the form given in Eqs. (33), (34) in the main text. 

It remains to interpret the right hand side (W ext ) of (A7). First, using explicit expressions (A2), (A5) for j, (j), 
we obtain (cf. (31)) 

(Wext) = J dS (wir 12 ) (A9) 

So 



{W ext ) = J dS\{-iuwf>6T$ +iuju { ° ) 8T[° ) ) (A10) 

So 

Eq. (A10) may serve to check equation (A7) explicitly. Two points appear explicitly : (i) Only those irradiated parts 
of the interface where iii ^ do contribute: the effect is connected with the free sliding of the interface points along 
the sliding direction, (ii) All three waves superimposed enter (W ex t). This is particularly remarkable in the case of 
an evanescent L wave, in which there is no longitudinal flow at infinity, yet the energy gain at the interface cannot be 
obtained correctly without including the L wave contribution right at the interface. 
At each point of the interface, the external force acting on medium M is 

F = -rn = T*n + Rc{5t)u = F*+5F 

The time average is 

(F) = F* 

The averaged macroscopic force acting on S is thus J~s = SF* . Medium M' is pulled at the overall velocity v, thus 
the power spent by this force is J~s-v, independently of the presence of oscillatory acoustic fields. It is now easy to 
see that 
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T s -v = J dS ((« - u)-F) + (W ext ) (All) 
s 

The first term is the power dissipated against the friction forces: v — ii is the local relative interfacial velocity, while 
F has only a frictional component along the interface. Equation (All) thus expresses the partitioning of the total 
work done per unit time by the external force into the (irreversibly) dissipated power and the net acoustic gain. 
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FIG. 1. Schematic representation of the system: a transverse incident wave impinges at incidence 9 onto the sliding interface, 
giving rise to two reflected and two transmitted waves. 

FIG. 2. Relative powers wij(6) (I, J — T,L) reflected from the interface between a compliant medium and a rigid one versus 
incidence angle 9. Friction is velocity independent. R = oo, / = 0.2, ct/cl = 0.5. Incident wave (a) transverse, wtt ( — ), 
wtl ( ) (b)longitudinal , wll ( — ), wlt ( )• 

FIG. 3. Evolution of the relative reflected power w T t{9) with elastic contrast ratio R for Coulombic friction. 

/ = 0.2, c T /c L = 0.5, R — oo (• • ■•), 40 ( ), 20 (- ■ - •), 10 (— )• (a) Full incidence range; (b) Blow-up for 9 close to 

— 9lim = —30°. 



FIG. 4. Evolution of the relative reflected power wtt(9) with strengh A of the velocity depence of the friction coefficient 

for rigid M' . f = 0.2, c T /c L = 0.5, R = oo. A = (— ), 1 ( ), 2 (- • - -).(a) Full incidence range; (b) Blow-up for 9 close 

to —9um = —30°. Two more curves, for A — 0.5 (■ ■ ■•), 0.9 (thin full line) show the shift and narrowing of the singularity as 
A — A cr — 1 is approached. For A > A cr , the peak amplitude becomes finite. 



FIG. 5. Evolution of relative reflected powers with increasing A. Incident T wave, R = 20, / = 0.2, ct/cl = 0.5. A = 
( — ), 1 ( ), 2 (— ■ - ■). (a) wtt{9), (b) wtl(9). The height of the peak for 9 = —9i im is maximum, but finite, for A — 1. 



FIG. 6. Same plot as in Fig. 5, but for an incident L wave, (a) wll{9), (b) wlt{9). 



FIG. 7. Static interface: relative powers reflected from incident wave (a) transverse, wtt ( — ), wtl ( ) (b)longitudinal , 

wll ( — ), wlt ( )• R= 20, c T /c L = 0.5 

FIG. 8. lastically symmetric case (R = 1), and Coulombic friction (^4 = 0): relative powers reflected from incident wave (a) 
transverse, wtt ( — ), wtl ( ) (b)longitudinal , wll ( — ), wlt ( )• / = 0.2, ct/cl = 0.5. 

FIG. 9. Evolution of the relative reflected power wtt{9) with strengh A of the velocity dependence of the friction coefficient 
for symmetric case R = 1. / = 0.2, c T /c L = 0.5. A = (— ), 1 ( ), 2 (■ - ■)• 



FIG. 10. Energy flows for rigid M' . 
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